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ABSTRACT 

This paper presents the results of the 2009-2010 monitoring sessions of the starburst 
galaxy M82, obtained with the Multi-Element Radio-Linked Interferometer Network 
(MERLIN) at 5 GHz and e-MERLIN at 6 GHz. Combining several 5 GHz MERLIN 
epochs to form a map with 33.0 /iJy/bm noise level, 52 discrete sources, mostly super- 
nova remnants and Hn regions, are identified. These include three objects which were 
not detected in the 2002 5 GHz MERLIN monitoring session: supernova SN2008iz, the 
transient source 43.78+59.3, and a new supernova remnant shell. Flux density varia- 
tions, in the long (1981 to 2010), medium (2002 to 2010) and short (2009 to 2010) term, 
are investigated. We find that flux densities of SNRs in M82 stay constant in most 
of the sample (~95%). In addition, aside from SN2008iz and the well-known variable 
source 41.95+57.5, two sources display short and medium term variations over the 
period 2009-2010. These sources being among the most compact SNR in M82, these 
flux density variations could be due to changes in the circumstellar and interstellar 
medium in which the shocks travel. 
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INTRODUCTION 



M82 is one of the closest (d = 3.2 Mpc; iBurbidge et al.l 
1964) starburst galaxies known, producing a large popu- 
lation of massive, rapidly evolving stars, and an equally 
large number of supernovae. The central regions of the 
majority of starbursts suffer from high visual extinction. 
Radio observations, which are unaffected by dust, are thus 
ideal for the detection and study of these supernovae or 
supernova remnants (SNRs). 

Over the past 30 years, M82 has been subject to fre- 
quent radio monitoring at cent imetre wavelengths w i th the 
Very Large Array (VLA; e.g. iKronberg et all l200d ; lAllenl 
1 19991 ) and the Multi-Element Radio-Linked Interferometer 
Netw ork (MERLIN; e.g. iMuxlow et al.|[l994l ; iFenech et~all 
2008). With the detection of over 50 discrete objects, in- 
cluding over 30 SNRs, these regular observation programmes 
have the advantage of tracking the evolution of supernova 
remnants as their shells expand, which provides a way to 
investigate properties of the Inter-Stellar Medium (see the 
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review bv lWeiler et~afll2002l h Since sources in M82 can be 
effectively treated as being at the same distance (their rela- 
tive distances differing by only a few hundred parsecs, which 
is small compared with the overall distance of the galaxy), 
they can be studied with essentially the same linear resolu- 
tion and luminosity sensitivity. Another advantage of study- 
ing extragalactic SNRs is that younger objects can be ob- 
served in M82 (from a few to a few hundred years) than in 
the Milky Way (until recently, the y oungest was Cassi opeia 
A, which is over 300-yr old; however. iGreen et al.l (|2008T ) dis- 
covered a younger galactic SNR ^150-yr old). In add ition, 
new sources such as SN2008iz (iBrunthaler et all I2009T ) or 
the radio transient 43.78+59.3 |Muxlow et a.l.ll2010h would 
not have been detected without regular radio monitoring. 
Finally, the number of SNRs detected provides a way of de- 
termining star formation rates in starburst galaxies, based 
on the number of sources in the galaxy, which is related to 
the supe rnova rate, and by assuming a give n initial-mass 
function jCondonll 19921 ; iBotticella et ai]|2012l ). 

Regular observation programmes also allow for the 
monitoring of flux density variability in sources (e.g. 
IKronberg etafl l200Ch . such as 41.95+57.5, which has 
shown a continued decrease in flux density since its first 
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Figure 1. Position of the 52 sources in M82 identified from the 5 GHz MERLI N combined data pr esented in this work (contour plot at 
-1, 1, 1.414, 2, 2.828, 4, 5.656, 8, 11.282, 16 X 99 fiiy/hm). Sources present in lFenech et al.l {2008) are represented by circles while new 
sources are represented by triangles. 
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Figure 2. Measured flux densities of the phase calibrator 
0955+697 over the period 2002-2010. 



Table 1. Details of 2009-2010 MERLIN (EP1-7), 2010 e- 
MERLIN (EP8) and archival MERLIN observations. ToS corre- 
spond to time on source (in hours) a and the rms noise is given in 
/ijy/bm 6 . Unless indicated otherwise, the observation were done 
without the Lovell telescope. 
T Tabley Telescope not available. 
D Dcfford Telescope not available. 
L Including Lovell telescope. 

The table also displays, in the last column, measured values of 
flux density (in mjy) c of the phase calibrator 0955+697. The 
flux density for EPS, measured at 6.0 GH z, was converted to 
4.995 GHz using S tx v a where o=-0.23 JVollmer et al.ll2010T) 
for this source. 



observation in 1965 (|Trotmanlll996l ; iMuxlow et al.ll2005h . 

This paper presents and discusses the results of the 
2009-2010 monitoring sessions, which includes seven epochs 
at 5 GHz with MERLIN and one epoch at 6 GHz with e- 
MERLIN. These results have been supplemented with re- 
sults from the literature (|Kronberg et al . 2000; Fcnccl Tet al.l 
2008) increasing the time-span of monitoring observations to 
30 years. It is expected that most sources will show no (or 
only very little) f lux density variations over th e time periods 
prese nted here (jUlvestad fc Antonuccil 1 19941 ; IWeiler et al.l 
|2002| ). This work also intends to investigate the behaviour 
of SNR in M82 and comp are it to previous models (e.g. 
ISeaquist fc Stankovicll2007r ). 

Details of the data reduction techniques are given in 
Section[2] while Section[3]describes the source identification 
process and flux density measurements. Finally, flux density 
variations in the identified sources are discussed in Section 
[4] including more detailed discussions about a few sources 
in Section |+T] 
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2 MERLIN OBSERVATIONS AND IMAGE 
PROCESSING 

The 2009 monitorin g campaign of M 82 consisted of seven 
wide-field MERLIN jThomas son 1986) observations at a fre- 
quency of 4.994 GHz, observed between May 2009 and April 
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2010. These were made using parallel hands of circular po- 
larisation, and were correlated with a total bandwidth of 
16 MHz divided into 32 channels. Across all epochs com- 
bined a total on-source integration time of 286.5 hr was used. 
Each observing epoch was reduced and analysed individu- 
ally and a deep exposure map was produced by combining 
all of these data. In addition to these, a single e-MERLIN 
observation was included in this study. This observation was 
made in December 2010 as part of the e-MERLIN commis- 
sioning programme and used a total bandwidth of 512 MHz 
with a median frequency of 6.26 GHz. These data were cor- 
related into four individual sub-bands each divided into 512 
frequency channels. These e-MERLIN observations were ob- 
served prior the installation of the new e-MERLIN wide- 
band IF system in spring 2011 and consequently only one 
hand of polarisation was used and the data displayed re- 
duced sensitivity in parts of the observing band. Further 
details of both the MERLIN and e-MERLIN observations, 
including archival data, are provided in Table [1] 

For both MERLIN and e-MERLIN data, the calibrator 
3C286 was used to set the primary flux density scales and 
the point source calibrator OQ208 was used to determine the 
bandpass and relative gains of the antennas. Observations 
of the phase reference source, 0955+697, were interleaved 
throughout each of the observations and used to determine 
the telescope phases. Following phase calibration these data 
were re-weighted appropriately to account for the relative 
sensitivity of each MERLIN antenna. All epochs, as well as 
the combined 2009 epoch, were then imaged separately. 



3 SOURCE MONITORING 
3.1 Source identification 
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iFenech et ail (120081 ) identified 55 sources during the 2002 
deep-integration MERLIN observation of M82, and 49 
of them, with Speak 99/iJy/bm (~ 3a), were de- 
tected in our 2009-2010 combined field. Of the other six, 
three sources lie outside of the imaged area (37.53+53.2, 
38.76+53.5 and 46.56+73.8) and three have S peah < 
3a (42.62 +59.9, 42. 66+51.6 and 43.55+60.0). Supernova 
SN2008iz (B runthaler et ai1l2009T ) and the transient source 
43.78+59.3 (|Muxlow et al.ll2010l ) were also identified. 

In order to search for possible new sources, each of the 

51 identified sources were imaged and subtracted from the 
uv-data. The resulting dataset was then re-imaged, using the 
same parameters as the original field. Following inspection of 
this new image, one new source was identified (44.28+61.1, 
see £|4.1[) . bringing the total number of identified sources to 

52 (Figure Q] and Table [2]). 



3.2 Flux density measurements 

The flux densities of all compact sources in M82, as well as 
of the phase calibrator 0955+697, were measured for each of 
the eight epochs using the AIPS task tvstat, and any back- 
ground flux density surrounding the sources subtracted. In 
order to convert all flux densities in this work to 5 GHz, 
spect ral index information f or each source was retr ieved 
from iMcDonald et al.l (|2002T ). iBrunthaler et al.l (|2009l ) and 
iMuxlow et al. l ipoid ). When unavailable, standard values of 
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Figure 3. Top: Light curve for the bright stable sources 
(from top to bottom) 40.68+55.1, 43.31+59.2, 44.01+59.6 
and 45.17+61.2 used fo r internal calibration, using data over 
the periods 1981-1991 l|Kronberg et al.ll2000h and 2002-2010 
(present work). Bottom: Magnitude of the internal calibration 
corrections applied. 



a = 0.0 and a = —0.7 were used for Hn regions and SNR 
respectively. 

The overall flux density calibration scale, as derived 
from observations of the standard flux density calibrator 
source 3C 286 on an epoch by epoch basis, is accurate 
to few percent, however some small systematic variations 
in the flux density scale may remain between individual 
epochs. Since the phase reference calibrator shows signifi- 
cant flux density variations on the timescales of this moni- 
toring project (Figure [2}, an intern al cross-calibration pro - 
cedure following the methodology of iKronberg et all (|2000| ) 
was applied to the results to assess and correct for any sys- 
tematic differences in the flux density calibration of these 
data. Following the initial calibration and imaging of these 
target data, the flux densities of the bri ghtest and m o st his - 
torically stable SNRs identified by IKronberg et al.l (|2000h 
were measured to assess whether any systematic, inter-epoch 
trends were visible in their flux density. 

The bright sources 40.68+55.1, 43.31+59.2, 44.01+59.6 
and 45.17 +61.2 which showed a consistent trend (as 
defined by IKronberg et all 120001 ) were used to provide a 
re-normalisation of the flux density scale between individual 
epochs (Figure [3]). For each source, the differences between 
flux density measurement at a given epoch and the best-fit 
line to the data is computed. The resulting corrections for 
each epoch was computed by averaging these differences 
between each source, epoch by epoch. 
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Table 2. Combined epoch (averaged over all 2009 data) normalised (as described in ^33} 2009 peak and integrated flux densities, 
for all 52 detected objects. S ources previous l y iden t ified in the l iterat u re as SNRs o r Hll regions have been labe lled accordingly. These 
identifications are taken from lMuxlow et all jl994h . IWills et all l|l997h . lAllenl l|l999h and lMcDonald et al.l 42002T I. The Right Ascension 
and Declination arc offset from 09h55m00s and +69°40'00"(J2000), while the source name correspond to the equivalent offset in B1950 

coordinates . Peak flux densities we re measured with an error of ±0.0 1 mjy. 

*SN2008iz l|Brunthaler et alj|200sh and transient source 43.78+59.3 l|Muxlow et al]|201Ch . 
t Source absent from the 2002 data with no previous reference found. 
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3 4. MERLIN a nd eMERLIN light curves for some sources in M82. The 2002 and 2005 data were retrieved from lFenech et aTl 
and lAreol <2006l) respectively. Flux densities have been normalised as described in t|3. 21 
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continued 
4 RESULTS 

Since very rapid changes in flux density are not physically 
expected for the observed sources, and epoch-to-epoch vari- 
ations can be strongly dependent on calibration methods, 
the approach to measure variability used in this work does 
not examine these very short-term changes. With these 
considerations in mind, we find that most of the sampled 
sources in M82 show no significant changes in flux densities 
over all timescales. 



Figure [4] - 

Details of the averaged flux density measurements (over 
the seven 2009 MERL IN epochs) are giv en in Table [21 In 
addition to the 2002 |Fenech et al.l 120081 ') and 2009 data, 
flux density measurements (which were also normalised as 
described ab ove) from a 2005 MERLIN monitoring session 
(|Argol l2006h were also included. These observations were 
conducted at a frequency of 4.754 GHz in December 2005, 
with parallel hands of circular polarisation, a bandwidth of 
16 MHz divided into 32 channels, and a total on-source in- 
tegration time of eight hours. 



Data for each epoch can be found in Appendix [A] and 
a few select light curves are displayed in Figure [4] 



To improve our analysis o f variation, data take n over 
an 11-year period from 1981 l|Kronberg et al l l200Ch were 
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Figure 5. Short (over the 2009-2010 period, left), medium (over the 2002-2010 period, centre) and long (over the 1981-2010 period, 
right) term variability, as described in jj4] Names correspond to sources with rates of change outside a 2.5c limit of the distribution. 
SN2008iz and the variable source 41.95+57.5 lie well outside the range of this plot and are not represented. 



Table 3. Characteristics of the sources showing flux density variations (shown in large font in the table), aside from 41.95+57.5 and 
SN2 008iz, in M82. Short term decrease rates were derived from the 2009-2010 data while sizes, radial velocities and ages were retrieved 
from iFenech et al Medium term decrease rate s were computed over the period 2002-2010 and long term decrease rates over 

1981-2010, where the 1981-1991 data were retrieved from lKronberg et al.l i2000f ). The corresponding light curves are shown in Figure[6] 
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added to the light curves when available (21 out of 52 
so urces). A systematic sh ift of ~1.07 mjy between the VLA 
of IKronberg et al.l (2000) and MERLIN data was removed 
from the VLA data. This difference of the overall flux density 
of each sou rce is due to the higher beam size of the VLA ob- 
servation of lKronberg et all (|2000t) . which will be sensitive to 
more diffuse matter than the MERLIN data presented here. 
Changes in the flux density of the detected sources were 
studied in the short, medium and long term (over the 2009- 
2010, 2002-2010 and 1981-2010 periods respectively). A line 
S(t) (with t in months for short term and years for medium 
and long terms) was fitted through every light curve and 
the slope used to determine which sources displayed consis- 
tent flux density variations. The distribution of slopes was 
fitted with a Gaussian model (Figure [SJ). Outliers (with rate 
of change outside of 2.5a) indicate sources with noticeable 
variability. Note that SN2008iz (X s hort=~ 0.9 mjy/month) 
and the variable source 41.95+57.5 (X; on9 =— 3.1 mjy/year), 
which show extreme variation, are not represented in these 
plots. In addition, the Gaussian distributions in Figure [S] 
show an offset from zero, possibly showing a consistent de- 
crease in flux densities of all sources at all epochs, although 
within one sigma of the measurements. 

The left panel of Figure [5] shows flux variability in 
the short term (over the period 2009-2010). The sources 
43.18+58.2 and 44.01+59.6 show decreases in flux density. 
The central panel of Figure [5] shows flux density variabil- 
ity in the medium term (over the period 2002-2010). The 
sources 44.01+59. and 46.52+63.9 show steady changes in 
flux density. Finally, the right panel of Figure [5] shows flux 



variability in the long term (over the period 1981-2010). No 
evidence of systematic variation in flux density is seen in 
any of the 21 sources considered. All variability values are 
shown in Table |3j while the light curves for varying sources 
are displayed in Figure O 

We note that most of the sources showing short term 
flux density va riation are the most compact objects on VLBI 
scales in M82 (|Fenech et al.ll2008f ) after 41.95+57.5, whose 
flux-density is known to steadily decrease by ~8 . 8%/y r 
jTrotmanlll9"96l ; IKronberg et~aHl2000l : iMuxlow et all 120051 ). 
Each source showing possible flux density variations is dis- 
cussed in more detail below. 



4.1 Variable sources 

A total of six sources have shown flux density changes (i.e. 
with rate of change outside of 2.5<r of the distribution) 
in the short, medium or long term, including 41.95+57.5, 
SN2008iz and the transient source 43.78+59.3. The latter, 
being quite particular objects, will be discussed separately 
in Sections 14.1.11 14.1.21 and 14.1.31 Some characteristics of 
the other sources, including flux density rates of change and 
sizes, are given in Table [3] 

As previously stated, two of the three other varying 
sources (43.18+58.2 and 44.01+59.6) have small angular 
sizes, and are in fact some of the most compact SNR in 
M82 after 41.95+57.5. Based on expans ion velocity m e asure - 
me nts from VLB! and M ERLIN data, IFenech et aH (|2008h 
and lBeswick et al. (2006) determined that these sources are 
amongst the youngest in M82, with ages of 55 and 140 yrs 
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Figure 6. Left: Short, medium and long term light curves over the p eriod 2009-2010 (top), 2002-2010 (centre) and 1981-2010 (bottom, 
where the 1981-1991 data were retrieved from Kronbcrg ct al. 2000]) for sources showing noticeable decrease in flux density (shown in 
bold font on the plots), aside from 41.95+57.5 and SN2008iz (Tabic [3J. For long term variation, flux measurements were averaged over 
EP1-6 for 2009 and EP7-8 for 2010 (as described in Table H}. A line Sit) = S + Xt was fitted to each curve (dashed line), where X 
characterises the variation in mjy/month and mJy/yr for short and medium/long terms respectively. Right: Contour plot at —1, 1, 
1.414, 2, 2.828, 4, 5.656, 8, 11.282, 16, 22.564, 32, 45.1286 X 99 /iJy/bm of the source from the combined 2009-2010 MERLIN data. The 
beam size of 50x50 mas is shown in the bottom left-hand panel. 



(as of 2002). Consequently, their short and medium term 
brightness variations could be explained by changes in the 
circumstellar and interstellar medium through which these 
shocks travel. This may imply that these sources are at a 
different temporal stage in their evolution compared with 
the older more stable sources. 

Finally, the third source showing flux density variation, 
46.52+63.9, shows only medium term variability. From its 
light curves in Figure |S] the variation can be explained by an 
anomalously low value of the source flux density measured in 
2002. The source is then excluded from the varying sample. 



4.1.1 41.95+57.5 



This source has shown a continued decrease i n flux density of 
^8.8%/yr since its fi rst observation in 1965 l|Trotmanlll996l : 
iMuxlow et al.| [2005). This decay rate would imply that, at 
birth, the source would have had a flux density of ~30 Jy 
(assuming a free e xpansion and a source age of ~80 years; 
iFenech et al.ll200sl ). These two facts (continuous decay and 
high flux density at birth) suggest that 41.95+57.5 is likely 
to be an exotic event and several suggestions have been made 
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as to its nature, including the pos sibility that it may have 
been an off-axis gamma-ray burst (|Muxlow et al.ll2005l ). 

As seen in Figure [7] 41.95+57.5 still shows a steady 
decrease in flux density. However, its rate of change in flux 
density has decreased from 8.8%/yr to 6.4%/yr in 1981 to 

6.4% in 2010, when fitting between 1993-2010. 

VLBI observations of 41.95+57.5 by iBeswick et al.l 

(2006) in November 1998 and February 2001 shows the 
source to be expanding at a rate of ~2000 ±500 km s~ . 
However, further VLBI data from 2005 (|Fenechl 120071 ; 
iFenech et ail |2010h shows that the apparent radius of 
41.95+57.5 may have decreased between 2001 and 2005 
(30.4 mas in 2001 and 26.85 mas in 2005). 

The unusually luminous radio supernova SN1986J in 
NGC891 and 41.95+57.5 show some distinct similarities. 
Both sources display an asymmetric radio structure on 
milliarcsecond scales, modest (few lOOOkms -1 ) expansion 
velocities, relatively high initi al radio luminosities an d a 
power-law flux density decay (jBietenholz et af] |2010| ). In 
recent VLA and VLBI observations the spectral index of 
SN1986J has been seen to be evolving rapidly above 5 GHz, 
deviating from the flux density decline observed at lower 
frequencies. This has coincided with 8.4 GHz and 15 GHz 
global VLBI observations detecting an e merging, inver ted 
spectrum source in the centre of the SNR (Bartel 200^). It 
has been asserted th at this new source may be a plerion, or 
pulsar wind-nebula (|Bietenholz et al. 1 120041 ). 

Both the possible decrease in ap parent angular size 
l|Beswick et al.ll2006l ; lF^nech et al.ll2008T ) and the small devi- 
ation from the historical flux density decline of 41.95+57.5 
in these recent 5 GHz flux density measurements could be 
consistent with the potential appearance of a new higher 
frequency central component within the shell, or that 
41.95+57.5 comprises of two components: a central, con- 
stant flux density one with a thermal-like spectrum, plus a 
larger component with a synchrotron spectrum that is simul- 
taneously expanding and declining in flux density. In both 
of these scenarios, up until recent times, the synchrotron 
source has dominated the observed characteristics of the ob- 
ject. As such, this source may be undergoing an evolution 
analogous to that observed in SN1986J, albeit in a much 
older remnant. Unfortunately, the MERLIN and e-MERLIN 
observations presented here do not have adequate angular 
resolution to image in detail the structure of this source and, 
as yet, higher frequency (>1.6 GHz) VLBI observations have 
not confirmed this hypothesis. Further analysis of the vari- 
ation in flux densities of 41.95+57.5 and SN1986J, both in 
time and frequency, would be needed to assess the similarity 
between the sources. 



4.1.2 42.81+59.6 - SN2008iz 

The light curve of SN2008iz from t he presented MERLIN 
data and supplemented by data from Marchili et al. 12010). 
is shown in Figure [S] Following iMarchili et al.l ( 2010h . the 
flux density decline after reaching peak brightness can be 
modelled using an exponential decay, 



The 2009-2010 data display a flattening of the light- 
curve which does not seem to be properly fitted by an ex- 
ponential model. A power-law decay of the form 

S(t) oc (t - t ) a (2) 

where a = —1.15+0.05, shows a better fit to the flux- 
densi ty decrease, wi t h a X r«rf~ 20-8 while the exponential 
fit of IMarchili et all (|2010l ). including the 2009-2010 data, 
yields x? ed =148.7. 

Whilst further later time measurements are required to 
confirm these results, there is tentative evidence that the 
radio light curve of SN2008iz is showing a small reduction 
in the rate of its flux density decline. 



4-1.3 43.78+59.3 - Transient source 

This source was not pres ent in the 2002 and 2005 data 
but was detected in 2009 (iMuxlow et al.ll2010l) and in the 
MERLIN observation presented here. iMuxlow et al.l (|2010h 
suggested that this object might either be due to accretion 
around a massive collapsed object or be a fain t and unusual 
supernova (Jo seph. Maccarone fc F ender 20l3)- In the latter 
case, we would expect the flux density to steadily decrease 
over time, as the source would be too young to have reached 
the remnant phase. However, in the observations presented 
here, it shows a flux density at 5 GHz of ~0.90 mjy, which 
does not seem to vary significantly over the period May 2009 
- December 2010. The supernova hypothesis is thus unlikely. 



s(t) = Ki(t-t y 



r = K 2 (t-t ) 



(1) 



where to is th e explosion date, wh ich was fitted to be Febru- 
ary 18, 2008 (|Marchili et al.ll2010l ). 



4.2 44.28+61.1- A newly detected SNR shell 

This source, located at position 09h55m53.00s 
+69°40'47.24" (J2000), was not visible in the 2002 MER- 
LIN data. It has a peak flux density value of S pea fe=0.18 
mjy/bm and an integrated flux density of 0.40+0.36 mjy 
from the combined field. Being shell like in shape (Figure[U]), 
it is most likely a supernova remnant. 

By measuring the source extent, the age of the SNR can 
be inferred. The radius of the source was determined using 
the AlPS annular profiles task iring, giving 7^140+5 mas 
(Figure [10]). At the distance of M82 (3.2 Mpc), this gives 
an estimated radius of 2.1 pc. lFenech et al.l (|2008l ) measured 
expansion velocities of SNRs in M82 to be between 2200 
and 10500 km/s, with a mean velocity of v exp = 5650 km/s. 
Using the latter value, we can estimate the age of 44.28+61.1 
to be ~360yrs. 



5 DISCUSSION 

5.1 Possible observational effects 

Before discussing possible physical explanations for the 
long-term observed flux density decay of the sources in 
M82 (or lack thereof), it would be prudent to consider 
any observational explanations fo r any a pparent variation. 
The data from iKronberg et all (|2000h are from VLA 
A-configuration observations of M82 with a resolution of 
350 mas compared to the 50 mas resolution of the 2009-2010 
MERLIN observations. As such, virtually all of the sources 
discussed that were observed in both sets of observations 
will be unresolved in these VLA data and resolved in these 
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Figure 7. Light curve for the source 41.95+57.5 over the period 
1981-2010 . Data for each period were retri eved as follow: 1981 - 
1991 from lKronberg et al.l (hoOOl) , 2002 from lFenech et al.l fcOOSl) , 
2005 from lArgol j200d) . 2009-2010 from this work. Models for the 
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iMuxlow et al.ll2005h and our best fit 6.4%/yr decay (solid line) 
for the period 1993-2010. The values of \ 2 quoted here correspond 
to the reduced chi-square statistics for each fit. 
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Figure 9. Contour plot at -1, 1, 1.414, 2, 2.828, 4, 5.656, 8, 
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Figure 10. Azimuthally averaged radial intensity profile and 
best-fit model (solid line) of SNR 44.28+61.1 obtained from 
MPS annular profiles task IRING. The peak intensity occurs at 
r^l40±5 mas, taken to be the radius of the shell. 
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Figure 8. Light curv e for SN2008iz. 5 GH z data from single dish 
Urumqi observations jMarchili et alj |2010) and the MERLIN/ e- 
MERLIN data are shown as open squares and filled circles re- 
spectively. An exponential de cay model (Equ . IT]| wa s fitted to the 
curve using parameters from iMarchili et al.l j2010h - solid curve 
- as well as a power-law decay model - dotted curve - which dis- 
plays the best fit to the decay part of the curve. The values of 
X 2 quoted here correspond to the reduced chi-square statistics for 
each fit. 

MERLIN data. As the MERLIN array lacks shorter base- 
lines, it is comparatively insensitive to the larger spatial 
scales, potentially 'resolving out' some of the emission 
associated with a given SNR. It is however considered 
unlikely that this would explain the observed flux density 
decay for the varying sources, and such an effect would 
become more pronounced for the larger, fainter and more 
diffuse SNR. All of the larger, diffuse sources do not show 
any flux density variation. This is compounded by the fact 
that the decay is not only apparent between the VLA and 
MERLIN data, but continues within the MERLIN-only 



data. 

The final observational consideration is that we are 
essentially measuring structurally evolving sources with a 
fixed system i.e. with a fixed uv-distribution which there- 
fore becomes sensitive to slightly different structures as a 
source evolves. However it is unlikely that this would ac- 
count for the flux density decreases observed, particularly 
over the timescales considered. All of the varying sources are 
well-defined, in general compact, and bright, thus such an 
effect would be more observable for higher velocity sources 
or those that are already more diffuse. 

5.2 Physical explanation 

Of the sources observed in the data presented here, two 
show short and medium term variations in flux density: 
43.18+58.2 and 44.01+59.6. These sources are some of the 
most compact objects in M82 after 41.95+57.5, SN2008iz 
and the transient source 43.78+59.3, with mea sured sizes of 
0.8 a nd 1.9 pc in MERLIN 2002 observations (jFenech et all 
l2008h . 
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As discussed in iFenech et alj l|201d\ the SNRs in M82 
are expected to be the result of core-collapse supernova 
events with a massive, typically red supergiant (RSG), pro- 
genitor star. At the time of the supernova explosion, there 
will be a complex circumstellar environment consisting of 
a slow- moving dense RSG wind with a density profile of 
p oc r~ 2 followed by a wind-blown bubble of almost constant 
low density produced by the main-sequence low-density, 
high- velocity wind. Beyond this wind-blown bubble, there 
will exist potentially the remain s of a collapsed H II region 
and eventually the ambient ISM (|Weiler et al.ll20o3 . and ref- 
erences therein). 

It is generally accepted that following the evolution 
through the wind-blown bubble, it is likely that there will 
be a peak in flux density as the supernova ejecta travels 
from the low-density wind-blown bubble into the higher 
density ISM (particularly in the case of M82). This will 
primarily be caused by the production of a thin-dense shell 
at the edge of the wind-blown bubble. There will however, 
also be a subsequent decline in flux density as the ejecta 
traverses this shell and emerges into the ISM, traditionally 
considered the transitional phase from a supernova to a 
supernova remnant. It is possible that these transition 
phases could provide an explanation for the observed flux 
density changes. 

M82 SNRs appear to be a more luminous population 
than those in the Milky Way and the LMC, and typically ap- 
pear to settle into a constant flux evolution on shorter time 
scales and at smaller diameters than their Galactic coun- 
terparts. This is commonly attributed to the higher density 
environment within the starburst region in which they are 
situated. The high pressures associated with starburst activ- 
ity reduces the ability of mass loss from the progenitor star 
to expand into the ISM. This implies that the denser the en- 
vironment, the shorter the distance over which a shock has 
to travel to have fully traversed the circumstellar medium 
and sweep up sufficient material to enter Sedov evolution 
and become a well established SNR. More extreme starburst 
galaxies, such as Arp220, have shown a mixed population of 
compact sources covering SNe, SNRs and transition objects. 
They all typically show smaller sizes and higher luminosi- 
ties, consistent with an extrapolation of evolution within an 
even denser environment than that fo und in M82. 

A study of SNRs in Arp220 by IBateiat et all (|201lh 
determined that 3-5 of 18 sources, with sizes between 0.2 
and 0.4 pc, show variations over the period 2006-2008 (i.e. 
there is a change in flux >50% between the two epochs), cor- 
responding to short timescale changes in this work. These 
numbers are in accordance with our findings here, consider- 
ing that Arp220 i s a more activ e starb urst galaxy than M82. 
From figure 5 of IBateiat et alJ (|201 lh . the variable sources 
in Arp220 and the compact variable sources in M82 fall into 
the range of sizes which is compatible with an SN to SNR 
transition (0.2-0.3 pc in Arp220 and 0.7-2.0 pc in M82). This 
prog ression could explain t he ch anges in flux density. 

ISeaauist fc Stankovid (|2007T j explore an alternative ex- 
planation for the compact non-thermal sources observed in 
M82, specifically that they could be wind-driven bubbles 
(WDBs) associated with very young super star clusters as 
opposed to SNR. One of the primary motivations for pre- 
senting this hypothesis is the lack of flux density variability 



for the m ajority of the compact sources within M82 as re- 
ported bv Kronberg et all (|2000h . The work presented here 
has shown that two of the most compact sources (apart from 
SN2008iz and 41.95+57.5) show flux density variability in 
the short and/or medium term over the period 2002-2010, 
all of which have measured expansion velociti e s of ^ 2700 
km/s |Beswick et alj|200d : IFenech et all [20081 . l20ld ), sug- 
gesting ages of ~ few hundred years. This can be added to 
the evidence provided by measured expansion velocities for 
many o f these sources which c ontrad icts the WDB hypoth- 
esis of ISeaauist fc Stankovid (|2007l ). The possibility that 
these compact sources are SNR in the transitional phase 
from circumstellar to interstellar medium fits well with a 
similarly observed population of sources within Arp 220, 
taking into account the higher pressures involved in the 
more intense star-formation environment within this galaxy. 
This would also suggest that the more diffuse non-thermal 
sources within M82 should show little variation in flux den- 
sity as they would now be expanding into the roughly con- 
stant density ISM. This does not rule out the potential for 
some of the non-thermal compact sources to be wind-driven 
bubbles and not SNR. It is interesting to note, however, that 
the two sources showing flux density variability with mea- 
su red expansion velocities, pass the feasibility tests outlined 
in ISeaauist fc Stankovid l|2007h , clearly suggesting that the 
proposed tests do not adequately distinguish between the 
WDB and SNR models. 

Continued regular monitoring of these sources will give 
us more insights to their behaviour. It is now possible to 
carry this out with new sensitive instruments such as e- 
MERLIN and the EVLA, and we expect to get a more com- 
plete view of the evolution of SNRs in M82. This will allow 
the investigation of fainter and/or transient objects in star- 
forming galaxies such as M82. 



6 SUMMARY AND CONCLUSION 

We presented the results of the 2009-2010 monitoring ses- 
sions of the galaxy M82, obtained with MERLIN and e- 
MERLIN at 5 and 6 GHz respectively. The data include 
eight epochs between May 2009 and December 2010, and 
were com pared with d a ta fro m the April 2002 monitoring 
session of IFenech et all (120081 ). 

A total of 52 sources were detected, including three 
sources absent from the 2002 data: supernova SN2008iz, 
the transient source 43.78+59.3 and a new SNR shell, 
44.28+6 1.1. Although the firs t two were previousl y ob- 
served l|Brunthaler et all 120091 : iMuxlow etahl |20ld ), the 
latter was first detected here. 

Light curves for all sources have been compiled, in- 
cludi ng flux density me asurements from 2002 (|Fenech et al.l 
l200Sf) . 2005 (|Argoll2006h and these new 2009-2010 data from 
5 GHz MERLIN monitoring sessions. We have found that, 
in most of the sample, flux densities of SNRs and Hn regions 
in M82 stay rel atively constant, in agreem e nt with previous 
analy s is fe.g. [Kronberg fc Sramekl 19921 ; iKronberg et al.l 
12000 ; IWills et all 11998m . Aside from SN2008iz and the 



well-known variable source 41.95+57.5, two sources show 
signs of changes in flux density in the short and medium 
term. These source are among the most compact SNR in 
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M82, possibly showing that flux density variations are due 
to changes in the circumstellar and interstellar medium in 
which the shocks travel. They could thus be transitioning 
from supernovae t o supernovae remnan ts, in analogy with 
recent results from iBateiat et all (|201lh in Arp220. 
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APPENDIX A: FLUX DENSITY MEASURMENTS 

Table Al. Normalised flux-density (as described in 33~2l for the 2002, 2005 and 2009-2010, MERLIN/ e-MERLIN monitoring sessions 
of M82, Peak an d integ rated flux densities are given in mjy/beam and mjy respectively. Spectral index information were retrieved from 
IrvlcDonald et all j2002T ). When unavailable, standard values of a = —0.7 for SNR and a = 0.0 for Hn regions were assumed. The 2010 
6-GHz e-MERLIN data were converted to 5-GHz using the relation S v oc u a . 
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